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We show a possible way to implement the Grover algorithm in large nuclear spins 1/2 < / < 9/2 
in semiconductors. The Grover sequence is performed by means of multiphoton transitions that 
distribute the spin amplitude between the nuclear spin states. They are distinguishable due to 
the quadrupolar splitting, which makes the nuclear spin levels non-equidistant. We introduce a 
generalized rotating frame for an effective Hamiltonian that governs the non-perturbative time 
evolution of the nuclear spin states for arbitrary spin lengths /. The larger the quadrupolar splitting, 
the better the agreement between our approximative method using the generalized rotating frame 
and exact numerical calculations. 



PACS numbers: 76.60.-k, 42.65.-k, 03.67.-a 
I. INTRODUCTION 

Recent experimentsBi have shown that electron spins 
in GaAs can preserve their coherence for distances of 
more than 100 /j,m and for times of up to 130 ns. Long 
coherence lengths and times are the main ijeaukjsment 
for performing logic (ffl^mtioiisdp^-s^ and 
quantum computing^maBEiaBm because both 
research fields are interested in the complete control over 
the phase information of the spins and qubits, respec- 
tively. Among the various available information carriers, 
nuclear spins have one of the longest coherence times due 
to their weak interaction with the environment, which 
amounts to storing phase information for a long time in 
e.g. semiconductor structures. Either tkgy can be used 
as qubits or in the unary representation.^ Although im- 
plementations based on the latter are not scalable, they 
are more feasible, more robust, and can also be imple- 
mented in classical systems.li3 Coherent access to the 
nuclear spins in semiconductors is achieved by the all- 
optical NMR method shown in Refs. |l^ and which 
makes use of the hyperfine interaction between the elec- 
tron and nuclear spins and relies on the large coherence 
time of the electrons. A further possibility to access 
the nuclear spins coherently is conventional NMR us- 
ing coils.Eil However, compared to the all-optical method, 
coils do not provide the spatially selective manipulation 
of nuclear spins. Recently incoherent transfer of elec- 
tronic spin to nuclear spin has been demonstrated exper- 
imentally in semiconductor structures u&kig,the quantum 
Hall effect or ferromagnetic imprint ing.E3E3 

One of the most interesting quantum algorithms was 
introduced by Grover ,Ej who demonstrated that the par- 
allelism of unitary operations can speed up the search 
for a desired quantum state. Since the Grover algorithm 
needs only the superposition principle of quantum or clas- 
sical mechanics, the unary representation of an ensemble 



of single particles-can be used, such as the atomic levels of 
a beam of atomsEj or the. large spin of molecular magnets 
embedded in a crystal.ES In this paper we show that the 
NMR method presented in Ref. |l^ works for arbitrary 
spin lengths /. We compute the Grover algorithm on the 
single-spin states |m) of the large spin of nuclear spins in 
semiconductors. For this it is essential that a delocalized 
state lip) = Y^ly^^_j CLrri \m) with arbitrary amplitudes 
can be produced, which we show to be feasible with a sin- 
gle magnetic rf pulse. Since arbitrary amplitudes CLfYi are 
needed, we derive a non-perturbative method to calculate 
the time evolution of the nuclear spins that is valid for ar- 
bitrary spin lengths /. We also perform exact numerical 
calculations based on solving the Schrodinger equation. 
It turns out that the larger the quadrupolar splitting, 
the better the agreement between our non-perturbative 
method and the exact numerical calculations. In order 
to be able to control the nuclear spins coherently, the 
energy levels Em have to be non-equidistant, which is en- 
sured by the quadrupolar splitting. Instead of encoding 
the information into the phases of a^, we encode the in- 
formation into the eigenenergies Sj-m of the eiffleijpt|ates 
|m) in the generalized rotating frame (see below) .EZIEj Af- 
ter choosing a specific basis state \M) to be looked for, 
we make it degenerate with the completely delocalized 
state \s) = cim 1^) that has equal amplitudes am 

in the generalized rotating frame, i.e. a second magnetic 
rf pulse lets \s) evolve into |M), which is ensured by the 
finite overlap of \s) and |M). 

In Sec. H we give the Hamiltonian that is the most 
suitable for nuclear spins in semiconductors. For the 
time evolution of the nuclear spins we first identify the 
transition (Feynman) diagrams that give the largest con- 
tribution to the multiphoton transition amplitudes by 
means of high-order time-dependent perturbation the- 
ory, which is demonstrated in Sec. HI. From this we 



derive effective Hamiltonians in the so-called generalized 
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rotating frame that describe the time evolution of the nu- 
clear spins non-perturbatively for arbitrary spin lengths 
/, which is shown in Sec. 0. In order to understand 
better our concept of the generalized rotating frame in 
Sec. 0, we first review the transformation to the standard 
rotating frame in Sec. In Sec. |v| we explain our pro- 
posed method to implement the Grover algorithm into 
the single-spin system of nuclei in semiconductors. Fi- 
nally, Sec. VII shows that conventional incoherent NMR 
can be used to read out the searched state after perform- 
ing the Grover quantum search algorithm. 
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FIG. 1. Multiphoton transition scheme (QC scheme) for 
the coherent population of the Iz eigenstates \m) of a nuclear 
spin 7 = 3/2. In this scheme the Grover algorithm can be 
implemented non-perturbatively. The frequencies Uk of the 
fields Hk are red (- •) and blue (- -) detuned. Diagrams con- 
taining blue detunings are negligible for large quadrupolar 
splitting, i.e. A ^ hSk > 0. In Sec. Ill we explain why only 
the leftmost diagram is relevant. 



II. MODEL HAMILTONIAN 



method shown in Ref. |2^ yields the following quadrupolar 
constants for Ga and As nuclei in GaAs semiconductors: 
A = 7x 10-^ K for ^^Ga, A = 3x lO'^ K for ^^Ga, and 
A = 2 X 10~^ K for '^^As. Let |m) be the eigenstates of 
7io with eigenenergies em- 
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2. Preparation of \s) = (1/^) Z)m=-i/2 1^) 
means of Eq. (^) in the QC scheme, which takes about 0.2 
ms for ifi = = 1 G, TTs = 0, = 6083 s'^ and 62 = 0. 
The analytical result is confirmed by numerics in Fig. I The 
calculation was done in the generalized rotating frame as de- 
scribed in Sec. 

Next, we apply external magnetic fields Hx^k{t) = 
Hk{t) cos{u;kt -\- ^k)^ k = 1,2,3,..., that differ by the 
phases and oscillate at frequencies uok being detuned 
by Sk from the differences between the energy levels 
Sm- To be more precise, cok is detuned by from 
£2/-/c+i — ^21 -k- Fig. |l| shows the example for / = 3/2 
and /c = 1,2,3. If(5/e<0 the detuning is said to be 
blue, and if > the detuning is said to be red. For 
GaAs, oOk-fOo'j^ ^ 10 MHz with 5k ^ ^ kHz, and a longitu- 
dinal magnetic field i^^ ~ 1 T is appropriate. It is desir- 
able to make Hz sufficiently large to accommodate many 
spin precessions before the spins dephase. The complete 
Hamiltonian has then the form 



The ensemble of nuclear spins in semiconductors is best 
described by the single-spin Hamiltonian 

'Hq = Hz + Hq^ (1) 

which consists of the nuclear Zeeman term 

Hz = —QnI^nHzIz, (2) 

with the nuclear abactor qn = 1.3 (see Ref. |l^), and the 
quadrupolar termEll 

Hq = A[3/2-/(J + 1)]. (3) 

The quadrupolar constant A differs significantly between 
the various nuclei. For example the all-optical NMR 



H = Ho^V{t), (4) 

where 

21 

V{t) = ^gN/J^NHk{t) cos{u;kt ^ ^k)Ix (5) 

k=l 

is the driving Hamiltonian given by the external magnetic 
fields. As usual, the spin operator can be decomposed 
into ladder operators, i.e. Ix = (/+ + /_)/2. 

For the implementation of |th|e-|non-perturbative ver- 
sion of the Grover algorithmEZlEj we have to be able 
to produce a completely delocalized state of the form 
\^) — Xlm=-/^^l^) with arbitrary amplitudes a^. 
Furthermore, we need to have coherent control over the 
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unitary time evolution of for arbitrary times t. We 

will show in the next sections that our non-perturbative 
method gains complete control over all the amplitudes 
am{t)^ which can be used in a future experimental imple- 
mentation. 
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FIG. 3. Preparation of \s) = (1/^3) El^=-i/2 1^) solv- 
ing the Schrodinger equation exactly for ^^Ga nuclei in the 
QC scheme, which takes about 0.2 ms for Hi = H2 = 1 G, 
H3 = 0, ^1 = 6083 s~^, and 62 = 0. The small oscillations 
are due to the five diagrams with blue detunings (see Fig. |^) 
that have been neglected in Fig. |. The ^^Ga nuclei have an 
anisotropy constant ofA = 3xlO~^K. 



III. PERTURBATIVE EXPANSION 

Usually perturbation theory is useful only if the trans- 
verse term V{t) is small compared to the longitudinal 
term Ho- But here we use high-order time-dependent 
perturbation theory to identify the transition diagrams 
with the largest contribution to the transition ampli- 
tudes Sm,m' = limt^_oo,t' ^00 {m\U{t,t')\m'). For this 
we first expand the ^-matrix in powers of V{t)^ i.e. 
^ — Xlj^o "^^^ jth-order term of the perturbation 

series of the ^-matrix in powers of the driving Hamilto- 
nian V{t) is expressed by 



dtk 



dtjQitk - t/e+i)/7o(oo, ti 



xV{h)Uo{tut2)V{t2) . . . V{tj)Uo{tj, -00), (6) 

which corresponds to the sum over all transition diagrams 
of order j, and where Uo{t, to) = e~*'^o*^^~^o)/^ is the free 
propagator, Q{t) is the Heavyside function. 

For illustration, we derive now a three a-photon tran- 
sition from |3/2) to |— 3/2). All the relative phases 
between the magnetic fields are assumed to be zero in 
this calculation. Then we obtain 



sn 3 = (ff ) 



dtidt2dt3e{ti - t2)e{t2 - ta) 



iesti/h —iei{ti—t2)/h —is_i{t2—t3)/h 

Xp_3 36262 62 
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i(a;_3 3—oJi—oj'2—oJ3)ti 



1 3 -u;i-u;2)(^i-t2) -iiuJi 3 -c^i)(^2-t3) /^n 
X6 2'2 e 2'2 (7j 

where from the 3^ = 27 terms remains only one due to 
the rotating wave approximation and after keeping only 
the terms with the smallest detuning. This diagram is 
the most left one in Fig. |l|. 




FIG. 4. Grover algorithm calculated by means of Eq. ( |65| ) 
in the QC scheme, where \s) — (l/y^) y^^^^_-i I'm) 
is concentrated mainly into | — 1/2) after 0.55 ms for 
H2 = n62/2gNl^N = 1 G, /ii = /i2, h3 = 0, 61 = 0. The 
duration of the QC is < l/2z^^^^^^-. This result is numerically 
confirmed in Fig. p^. The calculation was done in the general- 
ized rotating frame as described in Sec. 



Next, we make use of the relation 



■ i^m.mJ + + 



and substitute the Fourier transform 
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which yields 
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(10) 



as long as all Hk{^) have a sufficiently narrow peak 
around 1] = and Uk are detuned. 
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FIG. 5. Grover algorithm calculated by solving the 
Schrodinger equation exactly for ^^As nuclei in the QC 



scheme, where |s) = (1/ a/3) XI 



3/2 

m=-l/2 



\m) is concentrated 



mainly into | — 1/2) after 0.55 ms for H2 — ?i62/2gN fJ^N = 1 G, 
hi = h2, hs = 0, 61 — 0. The ^^As nuclei have an anisotropy 
constant ofA = 2xlO~^K, which is about 10 times larger 
than for ^^Ga nuclei. Therefore the small oscillations due to 
the neglected five diagrams in Fig. ^are here almost invisible. 
Hence we can infer that the larger the quadropular constant 
A, the better the agreement between our non-pert urbative 
method and exact numerics. 

After evaluation of the three-fold convolution we ob- 
tain 

q(3) ^ ^9}^Y J dtUk=i^k{t)e 



" 2 ' 2 \ ATi J ip \ s (^1 —^13) (ui -\- UJ2 — uj_i 3) 
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for 53 = 0, 
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for i52 = and H3 = 0, and 



(14) 



(15) 



for 5i = and H2 = Hs = 0, where 5* 
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FIG. 6. Grover algorithm calculated by means of Eq. ( |65| ) 

in the QC scheme, where |s) = {1/ V4:)^^J^'^_^^^\m) 
is concentrated mainly into |l/2) after 0.05 ms for 
hi — h2 — — h5i/2gN iJiN = 5 G, ^2 = = 0. 

It is interesting to note that 



lim ^^^^ — lim ^'^^^ 

^™ ^-3/2,3/2 - ™Q ^-1/2,3/2 



where Prn.m' = rife^m 1^- I A: + 1) is the spin amplitude 
between \k) and \k + 1) and 5^^\u) = ^ j^^j^ e'^'^dt = 
sin(a;T/2) ^j^^ dclta-functiou of width 1/T. In addition, 
we have used rectangular pulse shapes of duration T for 
all fields, i.e. 



0, (16) 
i.e. destructive interference is maximal. However, if 



A>n|4|, A: = 1,2,3, 



(17) 



likit) 



Hk for -T/2<t< T/2, 
otherwise. 



(12) 



The energy is conserved for uoT 1. Also, the duration 
T of the rf pulses must not exceed the dephasing time 
of the spin states. 

If we keep also the further five diagrams shown in 
Fig. [|, we obtain 



destructive interference is negligible. In other words, the 
most left diagram in Fig. |l| gives the largest contribution 
to the transition amplitude. This means that this dia- 
gram only is responsible for the time evolution of the nu- 
clear spins. So Eq. ( p!7| ) is the basis condition for finding 
the Hamiltonian in the generalized rotating frame that 
describes the time evolution of our nuclear spin system 
for arbitrary spin lengths / and arbitrary times t (see 
Sec.0). 
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IV. THE STANDARD ROTATING FRAME 

Before deriving the transformation of the Hamilto- 
nian TL including many magnetic fields Hx^kif) oscillat- 
ing at different frequencies uok to the generalized rotat- 
ing frame (see next section), it is instructive to have a 
look at the transformation of a Hamiltonian TYs includ- 
ing only a single circularly polarized oscillating magnetic 
field H^(t) = H±{cosu;tex + sin cotey) to the standard 
rotating frame. We start from the Hamiltonian 

= -QnI^nIHzIz + H±{Ix cos cot + ly sincot)] 

= -gN^HJ, ^ gNm\H^e'''''^he-'^''^ . (18) 

The proof for the second equality can be given as follows: 
First let us define the function 



The first and second derivative read 



-6 = -/(^)- 



(19) 

(20) 
(21) 



The solution of the last differential equation is 

f{(p) = Acosif ^ B sin (p. (22) 
The boundary conditions 

/(0) = A = 4, {^) ^__^ = B = -ly (23) 

reduce Eq. (H) to 

f{(p) = Ix cos (f- ly sin (f, (24) 
which proves Eq. (p!^). 



In order to derive the Hamiltonian in the rotating 
frame, we have to transform the Schrodinger equation 



ih 



d\ 



dt 



(25) 



with the unitary operation U. The left hand side of Eq. 
transforms into 



dU- 



dt V 9t dt 

and the right hand side of Eq. ( p5| ) into 



\Aot) 



(26) 



(27) 



where we have defined the state in the rotating frame 
by IV^rot) = ^l'^)- Combining both sides ( p6| ) and (|27| ) 
yields 



ih U 



dU- 



dt dt 
from which we obtain 

,^9|V^rot) _ f 



dt 



unu-^ - ihu^— 
V at 



iV^rot) = W-^lV^rot), (28) 



iV^rot) = Hrot iV^rot) 

(29) 



by multiplying U to Eq. (|2q ) from the left. From Eq. 
( p9[ ) we can immediately read off the Hamiltonian in the 
rotating frame: 



-Hrot = unu^ 



dt 



riTJ 

ih^Ul (30) 
dt ^ ^ 



If we insert U{t) = e 

Hrot = -{gNl^NHz - flUj)Iz + gNl^NH_iIx. (31) 

which is completely time-independent. In the next sec- 
tion we will show that under the special condition of 
quadratic anisotropy we can transform a Hamiltonian 
with more than one oscillating term into a Hamiltonian 
that is also completely time-independent. 
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FIG. 7. Grover algorithm calculated by means of Eq. ( pq ) 
in the QC scheme, where \s) = (1/V^) X^^=-3/2 1^) 
is concentrated mainly into |— 3/2) after 0.05 ms for 
hi — h2 — h2> — h52>/2gN iJiN = 5 G, = ^3 = 0. 
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FIG. 8. Grover algorithm calculated by means of Eq. ( |65| ) 

in the QC scheme, where |s) = (1/v^) Xlm=-5/2 1^) 
is concentrated mainly into | — 1/2) after 0.05 ms for 
hi — h2 — h^ — h4. — h^ — hS3/2gNlJ^N = 5 G, 
^1 = = ^4 = ^5 = 0. 
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V. THE GENERALIZED ROTATING FRAME 

For the case where the Hamiltonian Ti contains many 
magnetic fields Hx^k{t) that oscihate at different fre- 
quencies cj/c, the transformation to the standard rotating 
frame cannot be used anymore. In this section we develop 
a concept that allows us to transform the Hamiltonian 
H into a generalized rotating frame, i.e. H Hgrot, 
where TYgrot is completely time-independent. Then the 
time evolution of the nuclear spins can be calculated non- 
perturbatively in this generalized rotating frame. In con- 
trast to previous workEj our method also holds for van- 
ishing detuning energies hSk 0, which is essential to 
perform non-perturbative unitary operations. Once the 
control over 21 magnetic fields is established, the scheme 
proposed here allows for quantum information processing 
and quantum storage with a single pulse, provided that 
there is sufficient signal amplification due to the spin en- 
semble. The only requirements are that the quadrupolar 
splitting A is much larger than the detuning energies S^, 
so that only the most left diagram in Fig. |l| governs the 
time evolution of the nuclear spins (see Sec. |l|). 



7^(2) 



h2e-''^^^-'^^ £c 



(33) 



where hi = Hi {a |4| 6), h2 = H2 {b |4| c). 

The goal is now to find a unitary transformation U 

(2) 

that renders our Hamiltonian TYapp time-independent. 
We know from the previous section that the transformed 
Hamiltonian has the form 



(34) 

Instead of inserting U{t) = e"*'^^^^, we make an ansatz 
with three parameters cja, ^6, and coc' 



U 



-ioOat 










-iuJct 



(35) 




FIG. 9. Grover algorithm calculated by means of Eq 
in the QC scheme, where \s) = (1/^6)^ 



5/2 



-5/2 



is concentrated mainly into |— 3/2) after 0.05 ms for 
hi — h2 — — Ha. — — h54./2gN ijln = 5 G, 
5i — 82 — ^2, — 5b — 0. 

We first start with the example for two detuning en- 
ergies h5i and and give the general derivation for 
arbitrary spin lengths / afterwards. The corresponding 
Hamiltonian for three states reads 



7^(2) 



Sa {a\V\h) 
{h\V\a) Eb {h\V\c) 
(c|y|6) sc 



(32) 



where |a), and |c) are three consecutive spin states 
(a, 6, c are integers or half-integers with a = 6 + l = c + 2, 
— / < a^b^c < I) with their eigenenergies Ea^Sb^^c- For 
example a = 3/2, 6 = 1/2, and c = —1/2, which is shown 
in Fig. |l|. If the quadratic anisotropy A is much larger 
than the detuning energies hSk^ we can neglect all the 
transition diagrams with blue detuning energies shown in 
Fig. |l|. Applying also the rotating wave approximation 
leaves us with 
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FIG. 10. Grover algorithm calculated by means of Eq. ( |6q ) 
in the QC scheme, where |s) = (^/v^) Sm=-5/2 1^) 
is concentrated mainly into |— 5/2) after 0.05 ms for 
hi — h2 — hz — h4. — h^ — h5^/2gNiJiN = 5 G, 
5i — 82 — 5^ — 5a — 0. 

Evaluation of UH^pU^ yields 

Ea hie'^'^ 
hie-'^^^ eb h2e'^^^ 
h2e-'^^^ Sc 



(2) 
grot 



(36) 



with the abbreviations Qi = —uOa + cc;^ + cc^i + ^i/t and 
0.2 = —^b -\- ujc -\- ^2 -\- Since we want to make the 

Hamiltonian in Eq. (36) time-independent, we require 
that 



^1 = -OOa + Cc;^ + CJi + <I>i/t = 0, 
O2 = -^6 + C^c + ^2 + ^2/t 0. 



(37) 
(38) 
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In this way we get also rid of the Phases and ^2- 
Solving the Eqs. dH) and (H) gives 



^6 = ^2 + + 



(39) 
(40) 



ujc is a free parameter, which we choose to be ujc = —{uJi-\- 
iJ2)/2 - {^i + ^2)/2t. Then Eqs. (|3|) and ® simplify 
to 



(cJl+CJ2)/2+(^l+<^>2)/2t, 

(cj2-cJi)/2 + (<l>2-<^>i)/2t. 



Then the second term in Eq. (Q) is 



^^^^^ = 



fl(a;i+a;2) 
2 
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FIG. 11. Grover algorithm calculated by means of Eq. ( |6q ) 
in the QC scheme, where \s) = (1/vTO) y^^^_Q/2 1^) 
is concentrated mainly into | — 1/2) after 0.05 ms for 
hi — h2 — h2, — — hr, — Hq — h7 — hs — — MB/^gNfJ^N 
= 5 G, ^1 = ^2 = ^3 = ^4 = ^6 = ^7 = ^8 = ^9 = 0. 

If we insert the external frequencies uji = Sb/h — ea/h — 
Si and 002 = ^c/^ — Sb/fi -\- 5i — (^2, we obtain 



(2) 
grot 



- Sa-\-£c-h52 
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hi 




£a+£e+^(2(5l-(52) 
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h2 





/^2 
2 



which is equivalent to 



n 



(2) ^ 
^rot 



/ii 

/l2 ?i<52 



(44) 



(45) 



where we have subtracted ^^+^^ 



Id. Id is the identity 
matrix. 

Now we give the general derivation for the transfor- 
mation to the generalized rotating frame for arbitrary 
spin lengths /. From the above derivation, we know that 



for a spin of arbitrary length / the Hamiltonian in the 
generalized rotating frame has the form 



'^grot 



-me 



/lee'' 



/lee 



(46) 



which has non-zero entries only in the diagonal 
and first off-diagonal lines. |mo) = |2/) , |mi) = 
|2/ — 1) , . . . , |m2/) = |— 2/) are the eigenstates of 
the Hamiltonian shown in Eq. (|l|) with eigenenergies 
, ^mi , . . . , Sni2i • The Unitary transformation reads 



U 



with only diagonal elements that are non-zero. 



(47) 




FIG. 12. Grover algorithm calculated by means of Eq 
in the QC scheme, where \s) — (1/vTO) ^^^_g^2 1^) 



is concentrated mainly into 
hi — h2 — — h/i — — Hq - 
— b 5i — 82 — ^2, — — 5^ - 



1-3/2) after 0.05 ms for 
/17 = /i8 = /19 = M6/2gNfJ^N 
Sr = Sg = Sg = 0. 



For TYg^ot to be time-independent we have to solve 21 



linear equations of the form 



^2 

^6 

2/-1 



-cji + <l>i/t = 0, 
-002 + = 0, 



- ^5 

- ^6 



^5A = o, 

= 0, 



21-1 



21 



-^21 + ^21 /t = 0, 



(48) 
(49) 

(50) 
(51) 

A = 0, (52) 
(53) 
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Summing over all linear equations (Es^ to (p3|) leads to 



21 



, +^(Wfc + $fc/t), 



(54) 



k=l 



where u;m2i can be chosen arbitarily. Inserting this result 
back into the linear equations to (|53|) gives 



21 



^1712 — ^1712 



k=2 
21 



k=3 



21 



^1715 — ^1712 



— ^1712 



k=6 
21 



/c=7 



21 



, + ^ (wfe + <i>k/t). 



(55) 
(56) 

(57) 
(58) 

(59) 



k=2I 




UJQ 



jn- 



jn - ^6 + ^5 



^21 = £m2i/^ ~ ^m27-i/^ ~ ^21 + ^21-1 



(63) 



(64) 



After setting uor^ 



0, we then obtain for the diagonal 



elements of ih^U^ the following results: {—£mo +^ 



h6i - hd2i), {-e 



7712 



-772,27 



77727 ■ 



^77727 
dU 
dt 

^mi ~l~ ^77727 
^^2/)? • • • 7 ( ^m27-l ^ ^17121 

subtracting (£77727 ~ ^^2/)Id from TYg^ot Eq. ( [4 61) yields 
the time-independent Hamiltonian 



Sm^or + ^(^2/-i - ^^2/),0. Finally, 

(2/) 



(2/) 
grot 



hi 
hi hSi h2 
h2 hS2 



















h2i 
h2i h52i 



(65) 



This Hamiltonian allows us to evaluate the time evolu- 
tion of the nuclear spin system for arbitrary spin lengths 
/ and arbitrary times t non-perturbatively. Note that the 
Hamiltonian in Eq. (|65| ) remains valid even in the limit 
^ 0, where the perturbation expansions for multi- 
photon (more than 1) transition amplitudes, such as in 
Eqs. (^3|) and (p!^), break down. However, we must re- 
quire that \gN/-iNHk\ <C which means that the larger 
I A I, the faster the quantum information processing (see 
Ref. |l^). The time evolution of the state |'0(t)) reads 



(66) 



. -n / (27) , / 1- 

Propagators of the form [/T(t)e~ s^°* ' have 21 phases 
and 21 detunings h5k^ which determine the 21 phases 
and the 21 moduli of am- We subtracted two degrees of 
freedom: the global phase and the normalization condi- 
tion, respectively. In this way we can produce a state 
\ip{t)) with arbitrary amplitudes am- 



FIG. 13. Grover algorithm calculated by means of Eq. ( pq ) 
in the QC scheme, where \s) — (l/vTO) y^^^^_Q/^ \m) 
is concentrated mainly into |— 5/2) after 0.05 ms for 
/ii = /i2 = = ^4 = /i5 = /ie = /i7 = /is = /ig = ?i57 /2gN IJiN = 5 
G, ^1 = ^2 = ^3 = ^4 = ^5 = ^6 = ^8 = ^9 = 0. 



0.2 



Then the second term in Eq. 
the diagonal elements h{uOr) 



ih^W has 



diagonal elements of ih^U^ are zero. The external fre- 
quencies are 



All the off- 



/ 



\ 




- 197^1' 



ms 



^1 £mjh - Emo/h - Si 

^2 £m2/^ - £mjh -62^61 



^5 



5/^ - Sm^/h - + (^4 



(60) 
(61) 

(62) 



FIG. 14. Grover algorithm calculated by means of Eq. ([Gq ) 
in the QC scheme, where \s) = (l/V^) J^^=_9/2 1^) 
is concentrated mainly into |— 7/2) after 0.05 ms for 
hi — h2 — — Ha — — Hq — — hs — — h5s/2gNlJiN — 5 
G, ^1 = ^2 = ^3 = ^4 = ^5 = ^6 = ^7 = ^9 = 0. 
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VI. THE GROVER ALGORITHM 

Now we make use of the non-perturbative calculation 
method shown in the previous section to compute the 
Grover algorithm. We start from a configuration where 
mainly the ground state |3/2) is populated, sea. Fig. 
This can be achieved by the Overhauser effectEj. Then 
we produce an equal superposition of the nuclear spin 
states \s) = (1/v^) Zlm 1^)' which represents the ini- 
tial state of the Grover algorithm, n is the number of 
basis states involved in the search. For example we can 



3/2 



-1/2 



\m) with the 



prepare the state 1^) = (1/a/3)X1' 
parameters Hi = H2 = 1 G, i^s = 0, = 6083 s~^, 
and (^2=0, which is shown in Fig. ||. In Fig. || we cal- 
culated the time evolution of \s) exactly by solving the 
Schrodinger equation numerically. This calculation re- 
veals small oscillations that are mainly due to the five 
neglected diagrams in Fig. 



We continue the above example, where we have pre- 
pared the state \s) = (1/a/3) X]m=-i/2 1^)- ^^^^^ mak- 
ing \s) degenerate with | — 1/2) in the generalized rotat- 
ing frame, i.e. (^s H^l,^ 5^ = h52, the amplitude con- 
centrates mainly into the state | — 1/2) after 0.2 ms. In 
contrast to Refs. the Hamiltonian in Eq. ( |65|) has 

only nearest-neighbor couplings, which results in a de- 
creasing amplification of \M) with increasing / or |M|. 
However, even for the largest nuclear spin / = 9/2, we 
find that the resolution for identifying \M) is still suf- 
ficient, i.e. greater than 10% (see below). Again, we 
have also calculated this time evolution exactly by solv- 
ing the Schrodinger equation numerically, which is shown 
in Fig. ||. The small oscillations due to the five neglected 
diagrams in Fig. Q are now almost invisible. Thus the 
larger the quadrupolar splitting A, the better the agree- 
ment between our non-perturbative method and exact 
numerics. 




^/[ms] 



FIG. 15. Grover algorithm calculated by means of Eq. ( pq ) 
in the QC scheme, where \s) — (l/>/TO) ^^^_g^2 1^) 
is concentrated mainly into |— 9/2) after 0.2 ms for 
hi — h2 — — h/^ — — Hq — h7 — hs — — ?i5Q/2gN IJln - 
G, ^1 = = ^3 = ^4 = ^5 = ^6 = ^7 = ^8 = 0. 

Next, we choose one basis state \M) to be the one 
that we are looking for. Since we encode the information 
into the eigenenergies 5i-m of the eigenstates |m) in the 
generalized rotating frame, we can choose 5i-m = for 
all m, except for M, i.e. 5i-m 7^ 0. In order to find 
|M), we make it degenerate with the state \s) in the 



{21) 



grot 



^0. 



generalized rotating frame. So 5i-m - 

As we want to obtain the highest speed for the quantum 
information processing, the best choice for the Zeeman 
fields is to make them all equal, i.e. /ii = /12 = • • • = 
h2i =: h. The degeneracy condition now leads to 



{21) 
grot 



= -(2n- 
n 



2)h- 



1 2/ 

... 



k=l 



from which we obtain directly 

hSi-M 
h = . 



(67) 



(68) 



Energy 



-3/2) 



1-1/2) 



-|l/2) 



-|3/2) 



FIG. 16. In this example only the state |l/2) is fully popu- 
lated after performing the Grover algorithm. For the readout 
of this state, we have to irradiate the nuclear spin system 
with magnetic fields that oscillate at frequencies that match 
exactly the level separations. 

We have also calculated the Grover sequence for spin 
lengths up to / = 9/2. In all the cases for all spin 
lengths up to / = 9/2 and for all the searchable states 
—I<M<I the Grover sequence leads to an amplifi- 
cation of the amplitude of \M) between about 10% and 
100%. In Figs. I and the spin has length / = 3/2 
and the searched state is |M = 1/2) and \M = —3/2), 
respectively. When we increase the spin length /, the 
most interesting elements for applications in semicon- 
ductors are ^^Al, ^^Mn, and ^^Zn with nuclear spin 
/ = 5/2, and '^^Ge and ^^^In with nuclear spin / = 9/2. 
Therefore we have computed the Grover sequence for 
/ = 5/2 and M = -1/2, -3/2, -5/2, which can be seen 
in Figs. ^, ^ [1^, respectively, and also for / = 9/2 and 
M = -1/2,-3/2,-5/2,-7/2,-9/2, which is shown in 
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Figs. |ri|, |T2|, respectively. Because of the 

symmetry of the Hamiltonian in Eq. ^E^) along the di- 
agonal, the two computations for ±M look always the 
same. Hence, we have shown the Grover sequence only 
for half of the states -/ < M < 0. 



VII. READOUT OF THE RESULT 

After we have performed the Grover quantum search 
algorithm, we wish to extract the final state that we have 
been searching for. This can be achieved by using con- 
ventional pulsed NMR, where no coherence is required. 
Fig. [l6| demonstrates an example where only the state 
1 1/2) is completely populated. If we irradiate the sample 
with magnetic fields Hx^k{t) with zero detuning energies 
Sk = V/c, where the external frequencies are tuned in 
such a way that they match the energy level spaces, i.e. 
ujj-m = £m-\-i — ^m, then we induce only emission from 
1 1/2) to |3/2) and absorption from |l/2) to I -1/2). The 
emission and absorption spectrum shown in Fig. ^ iden- 
tifies unambigously the state |l/2). In this way we can 
read out any information stored in the nuclear spin sys- 
tem. 



using the effective Hamiltonian and exact numerical cal- 
culations. So the larger the symmetry breaking of the 
energy spacings, the better the control over the coherent 
time evolution of the nuclear spin system. 

The first test for the experimental feasibility of our 
proposal would be the implementation of multiphoton 
Rabi oscillations, as proposed in Ref. [l^. This can be 
achieved by applying only one oscillating magnetic field. 
Fig. ^ shows a 2-photon Rabi oscillation between the 
statesj3/2) and | — 1/2) for a nuclear spin 7 = 3/2. In 
general, multiphoton Rabi oscillations can be thought of 
as nutation of the large spin between the spin states |m). 

Once our scheme works experimentally, it would be in- 
teresting to apply a longitudinal magnetic field that has 
a large gradient. In this case the semiconductor sample 
could be divided into several regions, each of which could 
be addressed separately by different transverse magnetic 
fields. In this way one could perform parallel quantum 
information processes in each of the regions. Instead of 
using a large gradient field, it would also be possible to 
vary the g-factor in the sample, as was shown in the ex- 
periment of Ref. 30 for the electron spin in A^Gai-^^As, 



where the variation of the concentration of Al x leads to 




fA3/2, 1/2 CO-i/2, 1/2 



FIG. 17. The absorption and emission spectrui 
then unambigously the state |l/2) for the examp 
Fig. 0. 

VIII. CONCLUSION 

In our theoretical proposal we have showi 
Grover algorithm can be performed in the ni 
system of semiconductors. The main requirem 
Grover algorithm to work is the anisotropy i 
system, which is provided by the quadrupola 
The reason is that the quadrupolar splitting 
energy levels non-equidistant, which renders t 
spin states distinguishable. Only then we can find an ef- 
fective Hamiltonian of the form shown in Eq. (|65| ) that 
describes the coherent time evolution of the nuclear spin 
ensemble for arbitrary spin lengths and times. It turned 
out that the larger the quadrupolar splitting, the bet- 
ter the agreement between our non-perturbative method 



tSz 




FIG. 18. 2-photon Rabi oscillation between the states |3/2) 
and | — 1/2) for a nuclear spin 7 = 3/2. 
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